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We theoretically investigate the possibility of generating pulses in an excitable (asymmetric) 
semiconductor ring laser (SRL) using optical trigger pulses. We show that the phase difference 
between the injected field and the electric field inside the SRL determines the direction of the 
perturbation in phase space. Due to the folded shape of the excitability threshold, this has an 
important influence on the ability to cross it. A mechanism for exciting multiple consecutive pulses 
using a single trigger pulse (i.e. multi pulse excitability) is revealed. We furthermore investigate 
the possibility of using asymmetric SRLs in a coupled configuration, which is a first step toward an 
all-optical neural network using SRLs as building blocks. 
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I. INTRODUCTION 

Excitability is a phenomenon which is observed in a 
wide range of systems and has first been coined in bi- 
ology to describe the behavior of individual nerve cells. 
Excitable behavior is, however, by no means limited to 
nerve cells. Other examples include nonlinear chemical 
reactions, cardiovascular tissues, ion channels, climate 
dynamics and lasers pQ. Common to all these excitable 
systems is their highly nonlinear response to external per- 
turbations. When unperturbed, the system remains qui- 
escent and resides in a resting state. Small perturbations 
only lead to a small-amplitude linear response. How- 
ever, if the perturbation is sufficiently large, the system 
is transferred from the resting state to an excited state 
(the firing state). After this strong response, the system 
returns to its initial resting state through a refractory 
cycle. This large excursion of the system's variables in 
phase space corresponds to the emission of a large ampli- 
tude pulse. During the refractory cycle it is impossible 
to generate a second pulse — the system does not respond 
to any external perturbation. 

In optical systems, excitability has attracted much in- 
terest in recent years [2TfTT] . It provides a way to generate 
well-defined optical pulses and opens up possibilities for 
optical neural networks, all-optical pulse reshapers and 
delay lines. Lasers with saturable absorber [H [B], op- 
tically injected lasers [HI EH U and lasers with optical 
or opto-electronic feedback [2] [3] [5] [71 [9] have all been 
proposed as optical excitable units. 

In Ref. [H] we proposed a mechanism for excitability 
in systems with a weakly broken Z2-symmetry close to 
a Takcns-Bogdanov bifurcation [13]. As optical proto- 
types of such systems we used semiconductor ring lasers 
(SRLs), whose active cavity has a circular geometry. As a 
result SRLs can generate light in two opposite directions 
referred to as the clockwise (CW) and the counterclock- 
wise (CCW) mode. The convenient device properties of 
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SRLs allow this optical excitable unit to be highly inte- 
grable and scalable [HI [15] , allowing for fully integrated 
optical neural networks and all-optical devices. 

It was experimentally shown that short deterministic 
pulses can be excited by noise in asymmetric SRLs |12] 
and their origin was explained as a noise-activated escape 
across a barrier in an asymptotic 2D phase space [16] . 
Breaking the Z 2 -symmetry was realized by bringing a 
flat facet fiber very close to one of the output waveguides, 
increasing the reflective coupling for one of the counter- 
propagating modes. Moreover, that output waveguide 
was pumped to increase the power which is reflected [see 
[12] for more details]. Hence, an asymmetry in the mag- 
nitude of the linear coupling between the counterpropa- 
gating modes was obtained. 

From an application point of view it would be desir- 
able to excite pulses in a deterministic way by injection 
of an external optical trigger, which is theoretically inves- 
tigated in this article. We will also investigate the pos- 
sibility of forming an optical neural network using SRLs 
as building blocks. Optical neural networks are attrac- 
tive because of the high degree of parallelism that can be 
achieved and the large optical bandwidth that allows for 
very fast processing [T7] [TH] • The typical spike duration 
in asymmetric SRLs is of the order of 10 ns, which is 5 or- 
ders of magnitude faster than the biological ms timescale. 
Kravtsov et al. [18] have recently demonstrated an op- 
tical spiking neuron operating at GHz speed, but it re- 
quires a substantial setup. We feel that the ability of 
integrating the SRL based neuron on chip offers a clear 
advantage. These advantages in speed and size yield good 
perspectives as an artificial neural network. It will be 
shown in this article that coupled asymmetric SRLs are 
able to excite pulses in each other, mimicking neuron 
functionality as optical spiking neurons. 

This article is organized as follows. The modeling of 
the device — both the rate equation model and the asymp- 
totic phase plane — is discussed in Sec.[H] in which we also 
briefly review the excitability mechanism. In Sec. |III| we 
analyze the SRL response to optically injected pulses, 
revealing a sensitive dependence on the phase of the op- 
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tical trigger signal and a mechanism to excite consecutive 
pulses using a single trigger pulse. The possibility of cou- 
pling the SRLs toward the realization of an optical neural 
network will be examined in Sec. |IV| Finally, in Sec. [V] 
we summarize the results. 



II. MODEL 

We use a general rate equation model as proposed in 
[1"9] which assumes that the SRL operates in a single 
transverse and single longitudinal mode and can sustain 
two counterpropagating directional modes. It consists of 
two complex mean-field equations for the counterpropa- 
gating modes E\ (clockwise) and E% (counterclockwise) 
and a third equation for the carrier density N. Follow- 
ing [T5] with a straightforward modification to account 
for the optical injection, we can write the following rate 
equations for an optically injected asymmetric SRL; ne- 
glecting spatial variations within the laser and adiabat- 
ically eliminating the medium's polarization dynamics: 

£?i = k (1 + ia) [giN -1]E 1 - k 1 e i<t ' k E 2 + F(t) (la) 
E 2 = k (1 + ia) [g 2 N - 1] E 2 - h a t^"E 1 (lb) 



N = 7 



fi-N-g 1 N\E 1 \ 2 -g 2 N\E 2 \ 2 



(lc) 



Here the dot represents differentiation with respect to 
time t, E\ and E 2 are the slowly varying complex en- 
velopes of the counterpropagating waves, N is the car- 
rier population inversion, g\ = 1 — s \Ei\ 2 — c\E 2 \ and 
g 2 = 1 — s\E 2 \ — c \Ei\ are the differential gains of the 
modes (s and c respectively model the self- and cross- 
saturation effects), fj, is the renormalized injection cur- 
rent (/j, = at transparency and /i = 1 at lasing thresh- 
old), k is the field decay rate, 7 is the carrier decay rate 
and a is the linewidth enhancement factor. The term 
F(t) = Ei(t)e^ At+v ' /ri n represents the optically injected 
trigger pulses in the clockwise mode, where Tj n is the 
cavity round-trip time, Ei(t) 2 the power envelope of the 
injected pulse, A the detuning between the frequency 
of the injected field and the cavity resonance frequency 
and (p represents a constant phase difference. A detun- 
ing A > corresponds to an injected frequency which is 
higher than the cavity resonance frequency. The linear 
coupling between the counterpropagating waves, referred 
to as backscattering, is caused by reflections inside the 
cavity at the interface with the coupling waveguide and 
at the cleaved end facets of the output waveguide. It is 
modeled by a backscattering amplitude fej with a phase 
shift 0^. The backscattering phase <f>k is chosen to be 
identical for both counterpropagating modes since it has 
been shown that an asymmetry in has no influence on 
the topology of the phase space structure [TB] ■ The asym- 
metry in the backscattering amplitude Afc = k 2 — k\ is 
defined by the dimensionless parameter 8 — Ak/2k where 
k is the average backscattering amplitude. This asymme- 
try causes the SRL to be excitable [T^l [TB] . When 5 > 



in Eq. (JlJ, residence in the E 2 mode is favored, allowing 
for excitable pulses of the E\ mode. 

In a typical SRL, the photon lifetime k _1 and the car- 
rier lifetime 7 _1 are respectively of the orders 10 ps and 
5 ns, yielding two different time scales in the system. The 
other parameters are fixed to realistic values a = 3.5, 
s = 0.005, c = 0.01, the average backscattering value 
k = 0.44 ns -1 , 4>k = 1-5 and r, n = 0.6 ps [TS], un- 
less mentioned otherwise. The value of the bias current 
/i = 1.65 is chosen slightly above the value for which 
alternate oscillations disappear in a fold of cycles [T5] . 

Although we numerically investigate excitability in 
asymmetric SRLs with the rate equation model given in 
Eq. ([!]), the excitability mechanism can be interpreted 
more easily in a reduced two-dimensional phase space. It 
has been shown that on time scales slower than the re- 
laxation oscillations the dynamics of the SRL essentially 
take place in a two-dimensional phase plane. The dy- 
namical behavior in this phase plane is described by an 
asymptotically reduced model which has been introduced 
in 20], characterized by the variables 9 and ip. The re- 
sulting asymptotic description of the SRL is valid on time 
scales slower than the relaxation oscillations and it has 
been shown that it is able to predict many of the exper- 
imentally observed SRL characteristics [12l US HU [22] . 
But including optical injection in the asymptotic model 
makes it rather cumbersome and its validity can be ar- 
gued for short pulses. However, projecting our simulation 
results from Eq. (JlJ on the asymptotic phase plane will 
prove to be useful. 

The two phase space variables 9 E [— tt/2,tt/2] and 
ij) € [0, 2n] are defined by 



9 = 2 arctan 

1p = (j) 2 ~ <f>\. 



\Ei\ 



2' 



(2a) 
(2b) 



where <pi = arg(i?i). Hence 9 is a measure for the relative 
power distribution among the counterpropagating modes 
(6> = tt/2 if \Ei\ = 0, 9 = -tt/2 if \E 2 \ = and 6 = 
if \E\\ = \E 2 \) and ip is the relative phase difference be- 
tween the corresponding electric fields. The phase space 
topology of the asymmetric SRL in the (9, ip) phase plane 
is shown in Fig. [ljb). 

The grey and white regions indicate the basins of at- 
traction of the two stable states in the SRL — the CW 
and CCW state — which are quasi-unidirectional. In this 
case, the CCW state is favored. They are separated by 
the stable manifold of a saddle point indicated by S. The 
branches of the unstable manifold respectively spiral to- 
ward the CW and the CCW state. Assume that the un- 
perturbed SRL resides in the CCW state. A perturbation 
forcing it to cross both branches of the stable manifold 
will make the SRL relax back to the CCW state by turn- 
ing around the CW state, generating a large amplitude 
pulse. Hence the excitability threshold is defined by the 
stable manifold of a saddle point which is folded through- 
out the phase space and separates the basins of attraction 
of the counterpropagating states. 
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fashion. 

From the trajectory in the (8,xp) phase plane shown 
in Fig. [l|b), we see that the trigger pulse kicks the SRL 
across both branches of the stable manifold of the sad- 
dle (grey region). When the trigger pulse is finished at 
the point r the SRL relaxes back to the CCW state by 
turning around the CW state, which it is forced to do 
by the folded saddle manifold. The oscillatory tail of the 
pulse corresponds to the spiraling movement towards the 
CCW state in the (8, ip) plane and indicates the existence 
of subthreshold oscillations. While decaying toward the 
resting state (CCW), the close proximity of the stable 
manifold allows for some perturbations to experience a 
temporary reduction of the distance to the excitability 
threshold, which has been shown to lead to noise-excited 
double pulses P^lUB]. 



FIG. 1. (Color online) Simulation of Eq. |TJ when optically 
injecting a 2 ns wide square pulse, (a) Time trace of the 
modal intensities. The CW (CCW) modal power is depicted 
in red (black). The injected pulse is shown by a dashed black 
line. Note that the optically injected pulse power is scaled 
up by a factor 10 7 in this plot, r indicates the time at which 
the injected pulse ends, (b) Two-dimensional phase space 
trajectory corresponding to the time trace. The point r also 
corresponds to the moment when the injected pulse ends. S 
indicates the location of the saddle. The basin of attraction of 
the CW (CCW) state is depicted in gray (white). Parameter 
values: <5 = 0.045, Et — 8 x 10 -5 , resonant detuning, phase 
difference = 1.3-7T. 



We stress that the same results apply for every dynam- 
ical system with an underlying Z2-symmetry close to a 
Takens-Bogdanov point providing the same sequence of 
bifurcations [13] (for more information on the bifurcation 
scenario specific to SRLs, we refer to [T5]). 



III. OPTICAL TRIGGERING OF PULSES 

The asymmetric SRL will be excited to fire a pulse if a 
perturbation forces it to cross both branches of the stable 
saddle manifold by altering the phase difference ip and/or 
the relative power distribution 9 between the counter- 
propagating modes. The stable manifold can be crossed 
by introducing spontaneous emission noise as shown in 
P2]. In this article we investigate the possibility of de- 
terministically crossing the stable manifold by optically 
injecting a pulse into the SRL. An example of SRL "fir- 
ing" triggered by an optically injected pulse is shown in 
Fig. [ija). We have injected a 2 ns wide square pulse 
in the CW mode at resonant detuning (the optical fre- 
quency of the injected field and the SRL fields are iden- 
tical). The SRL responds by emitting a large amplitude 
pulse in the CW mode, while the CCW mode acts in 
antiphase showing a power drop. Subsequently the SRL 
relaxes back to its former state in a decaying oscillatory 



A. Phase dependency 

It is clear from Fig. [ljb) that the direction in which 
the SRL is kicked out of the CCW state has an influence 
on whether or not the stable manifold will be crossed, 
and hence whether or not the SRL will fire a pulse. 

We examine the initial direction in phase space by look- 
ing at the optical injection terms in Eq. ([I]). Prior to 
the injection, the SRL will reside in the resting state for 
which the right hand sides of Eq. are zero in absence of 
the F(t) term. Hence, at the start of the pulse the optical 
injection term is the only nonzero term. For the square 
pulse injected in the CW mode, F(t) = E i e i( - At+V '> /r m , 
and at the start of the pulse we get the following equa- 
tions for the amplitude \E\\ and the phase (pi of the field 
of the CW mode: 



d\E 



— = — cos(<^i - At Q - ip), 



df 



dt 



1 Ej 

\E\ \ T in 



sm{4>\ — Ato ™ ¥>)■ 



(3a) 
(3b) 



Here to indicates the time corresponding to the start of 
the injected pulse. Equations (|3| can easily be inter- 
preted by considering either constructive ((pi — Ato — 
ip = 0) or destructive ((pi — Ato — ip = ir) interference. 
Reformulating Eqs. (|3| for 8 and ip using Eqs. Q yields 



d9 _ 

~dt ~ ~\EL t 
dip 1 E, 



2\E 2 



Ei 



., cos(0i - Ai - <p), (4a) 

z + \E 2 r n, 



dt \Ei\Tii 



'2\ T in 

sin (0i — Ato — ip). 



(4b) 



Equations Q show that the angle under which the tra- 
jectory leaves the CCW rest state is exactly the phase dif- 
ference (p\ — Ato — f (defined as the clockwise angle with 
the negative 9 axis), which is an uncontrollable quantity 
in a practical setup. This phase difference will play an 
important role in whether or not the SRL will be excited 
to fire a pulse. An unfavorable phase difference leads 
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FIG. 2. (Color online) Influence of the phase difference be- 
tween the injected field and the SRL field on the response of 
the SRL. A pulse with a fixed width of 2 ns is injected in the 
CW direction with varying phase difference and amplitude 
Ei . Grey indicates excitation of a single pulse while black in- 
dicates excitation of two consecutive pulses. Inset: Magnitude 
of the response peak vs. amplitude of the optical excitation 
pulse for a fixed phase difference of 37r/2 (corresponding to 
the arrow). Note the clear threshold behavior. 



phase space (the stable manifold of the saddle point), it 
isn't accurate to speak of inhibitory and excitatory per- 
turbations. But the folding is such that there does exist 
a "wrong" direction for perturbations, along which no 
pulse will be excited [see Fig. [2] . The cumbersome phase 
dependency is an aspect that will arise in every optical 
excitable system, when triggered by an optically injected 
pulse. The experimental results of Wiinsche et al. in [TJ, 
covering an optically excitable multisection laser, indeed 
show a distinct jitter in the spikes although nominally 
identical pulses were injected, which is attributed to the 
phase difference between the fields. 

The inset in Fig. [2] shows the distinct threshold behav- 
ior for raising injection amplitude when the phase dif- 
ference is kept fixed. Below threshold, the response am- 
plitude increases proportional to the trigger amplitude. 
Above threshold, it jumps to a much larger value and 
is less sensitive to the trigger amplitude. It was shown 
in |23j that there exists an inverse correlation between 
the response pulse amplitude and the response pulse du- 
ration. The delay between the response and the trigger 
pulse ranges from approximately 3 to 10 ns depending 
on the phase difference, since this determines where the 
threshold curve is crossed and hence also the length of 
the phase space trajectory up to the pulse maximum. 



B. Multi pulse excitability 



to a perturbation in the wrong direction and hence fails 
to excite the SRL. Based on the topology of the saddle 
manifolds in the asymptotic (0,ij)) phase plane, we can 
predict which phase differences will have difficulties trig- 
gering a pulse. Given that Eqs. |2]) are only first order 
approximations with limited validity in time, and that 
the phase difference corresponds to the clockwise angle 
with the negative 9 axis in the (8, t/j) phase plane, a coarse 
estimate would be that phase differences between and 
7r/4 tend to be detrimental for triggering pulses. 

This is confirmed in Fig. [2j showing the influence of 
the phase difference on the firing of the SRL by numer- 
ical simulation of the rate equation model in Eqs. (JfJ. 
We inject a square optical pulse with a fixed width of 2 
ns at resonant detuning and monitor the response. In 
the diagram, both the amplitude Ei and the phase dif- 
ference between the injected field and the SRL fields are 
varied. For low injection amplitudes, the SRL only fires 
a pulse when the phase difference is close to — 7t/2. In the 
(9, ip) plane this corresponds to a downward kick [com- 
parable to Fig. |TJb)] , which is indeed the shortest way 
to cross the stable saddle manifold. For higher injection 
amplitudes, the phase condition is less stringent but nev- 
ertheless does not allow firing when the phase difference 
is close to 7r/4, which corresponds to our previous reason- 
ing. We can conclude that the phase difference between 
the fields determines the direction of the perturbation. 
Due to the folded shape of the excitability threshold in 



In the region of higher injection amplitude, the black 
area in Fig. [2] indicates the excitation of two consecutive 
pulses, an example of which is shown in Fig. [3} Note that 
the threshold for this double pulse generation is formed 
by the second fold of the stable manifold. By changing 
the values of the pump current J and the backscattering 
phase one can increase the number of folds of the sta- 
ble manifold [TH [TC] . This leads to identical scenarios in 
which more than two consecutive pulses can be excited 
using a single trigger pulse, i.e. multi pulse excitability. 
Multi pulse excitability was first reported for an opti- 
cally injected laser by Wieczorek et al. [5], explained 
by the vicinity of n-homoclinic bifurcations resulting in 
n response pulses. Other reported mechanisms are the 
period-doubling of a limit cycle on which a saddle-node 
bifurcation takes place, resulting in double pulses |10) . 
and slow-fast dynamics with folded slow manifold and a 
ramped parameter [24J. The mechanism reported here 
is characteristic for systems with a phase portrait sim- 
ilar to Fig. |ljb) , i.e. Z2-symmetric systems close to a 
Takens-Bogdanov bifurcation [13]. The multiple folding 
of the stable manifold, i.e. the threshold curve, around 
the resting state allows for a deterministic multi pulse re- 
turn trajectory given that the perturbation is sufficiently 
large [see Fig. [3^b)]. 

Note that this multi pulse generation differs from the 
noise-induced clustering we previously reported in [12] , 
In that case the SRL is excited a second time while it is 
relaxing back to the CCW state close to the stable mani- 
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FIG. 3. (Color online) Simulation of Eq. |TJ when optically 
injecting a 2 ns wide square pulse with a higher amplitude 
than in Fig. [l] (a) Time trace of the modal intensities, (b) 
Two-dimensional phase space trajectory corresponding to the 
time trace. Parameter values: S = 0.045, E, — 1.2 x 10 -4 , 
resonant detuning, phase difference = 1.3-7T. Conventions as 
in Fig.0 



fold, resulting in two very similar pulses. The determinis- 
tic return trajectory in Fig.|3^b) consists of a first, sharp, 
pulse followed by a wider pulse. The width of the excited 
pulses can vary depending on the strength of the initial 
perturbation, but the second pulse will always be wider 
than the first pulse. This range of widths is due to the 
characteristic flow in the reduced phase space [23]. The 
occurrence of this event in the noise-excited case [12] is 
very unlikely because of the the stochastic driving force 
and the invariant flow of the system. 



C. General considerations 

Since the excitability mechanism emerges from an 
asymptotic model which assumes the carrier density N 
to be close to its threshold value, the scenario described 
in this article must assume the amplitude of the trigger 
pulse to be small [i.e. O(10 -4 )]. In this way, the SRL 
is confined on the (9, tp) subplane of its phase space [see 
Sec. [n]. Injecting large pulses highly perturbs the value 
of N and causes the SRL to leave this plane. Using such 
a large amplitude trigger pulse, it is hence possible to ex- 
cite a pulse if the SRL happens to relax back to the the 
other side of the stable manifold on the (9, ip) plane. This 
excitability threshold is, however, obscured since the SRL 
dynamics can no longer be clearly described in the (9, tp) 
plane and the excited pulses have a distorted pulse shape 
due to large amplitude relaxation oscillations. The fact 
that the carrier density is not involved in the excitable 
excursion contrasts with the excitability mechanisms re- 



ported for other laser systems [H IfjHTU] . This is due to 
the ability of redistributing optical power between the 
counterpropagating modes, due to the backscattering k. 
In other laser systems one needs to generate additional 
optical power for the pulse for which the carrier reservoir 
needs to be addressed. 

Besides the amplitude of the trigger pulse, the detun- 
ing between the optical frequency of the trigger pulse and 
the SRL fields also has an influence. In the case of res- 
onant injection it is physically meaningful to talk about 
constructive or destructive interference between the in- 
jected field and the SRL field. However, if there is a 
frequency mismatch the interference will be constructive 
at one point, but destructive at another. This results in 
a periodic instead of constant driving force in the (9, tp) 
plane, which is unfavorable to cross the threshold (i.e. 
the stable manifold). Since the relevant physical quan- 
tity is the product At of the frequency detuning A with 
time t, as long as A < 0(1/t) (with r the pulse width) 
the driving force will not vary significantly. In that case 
the scenario sketched above remains valid. This is also 
what we observe numerically. We find numerically that 
there is a tolerance on the frequency mismatch of the or- 
der of a couple hundred MHz [cf. the pulse width is of the 
order of 2 ns], depending on the trigger pulse amplitude. 
There is of course a balance to be considered between 
the trigger pulse amplitude and the detuning. When the 
frequencies of the injected field and the SRL field are 
largely detuned one can always use a shorter pulse with 
a higher amplitude, for which the previous considerations 
(of large pulse amplitudes) apply. 



IV. TWO COUPLED EXCITABLE SRLS 

In this section we will investigate whether excitable 
asymmetric SRLs are able to excite each other, paving 
the way to an integrated optical neural network. A de- 
tailed study about the dynamical behavior of coupled 
SRLs lies outside the scope of this work and will be dealt 
with elsewhere. We limit ourselves to excitable behavior 
and consider two SRLs coupled by a single bus waveg- 
uide as shown in Fig. [4] In this type of coupling, the CW 
(CCW) mode of SRL a (b) is coupled in the CW (CCW) 
mode of SRL b (a) . We therefore extend the model with 
two coupling terms: 

Eu = k(1 + ia) [giaN a - 1] E u - k x e^«E 2a 

+ F(t) (5a) 

E 2a = k (1 + ia) [g 2a N a - 1] E 2a - k 2 e i ' t "°E la 

-k c e^E 2b (5b) 

E lb = k (1 + ia) [g lb N b - 1] E lb - k 1 j+*E 2b 

-kJ*'Eu (5c) 

E 2b = k (1 + ia) [g 2b N b - 1] E 2b - k 2 e^ k E lb (5d) 

The rate equation for the carrier density is given by 
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FIG. 4. (Color online) Schematic representation of the cou- 
pling scheme for excitable asymmetric SRLs considered in this 
article. The CW (CCW) mode of SRL a (6) is coupled in the 
CW (CCW) mode of SRL b (a). The notch to the ring cav- 
ity is merely added as a visual indication of the asymmetry 
of the cavity. The curved black arrow indicates the stable 
resting state, the dashed red arrow indicates the propagation 
direction of excitable pulses. The injected trigger pulse (F) 
is shown by a straight black arrow. 



Eq. ( lc ) for both SRLs, with respectively N — N a and 
N = Nb. The term F{t) again represents the external 
triggering and is only present in the equations for SRL a 
[see Fig. [1] . Based on the results of the previous section 
we will use a trigger pulse with a favorable phase differ- 
ence of — 7r/2 w.r.t. the fields in SRL a [see Fig. 111. The 
coupling is modeled by the coupling amplitude k c and 
the coupling phase <p c , which is the optical phase accu- 
mulated by the field when traveling from SRL a to SRL b 
(and vice versa). Note that Eqs. ^ are invariant under 
the transformation 



9c H> 

4>ib >-> 



I>c + 7T 

i« T tt/2 



(6a) 
(6b) 
(6c) 



where <j>ij = arg(£ , y ) and i = {1,2}. The second and 
third equation signify shifting the phase origin in each 
SRL over tt/2, but in the opposite direction. The only dif- 
ference with the untransformed system is that the value 
of the relative phase differences between the fields of the 
SRLs (without taking (j) c into account) will be shifted 
by 7T. However, the dynamical behavior of the global 
system will be identical. Therefore, we only need to con- 
sider <fi c 6 [0, 7r] to grasp the behavior for all values of <j) c . 
All parameter values are identical to the previous section, 
except for the asymmetry 8. We raise its value to 8% in 
this section to increase the stability of the resting state. 
The phase space of the solitary SRL is qualitatively the 
same as sketched in the previous section up to 8 = 8.4%, 
where the saddle and the CW state disappear in a fold 
bifurcation [see Ref. |16j for more detailed information]. 

For the coupling amplitude k c we have chosen k c < k, 
which means that the reflective coupling inside each SRL 
has to be larger than the coupling to the other SRL. 
This is hard to achieve when the backscattering is small. 
However, the value of the backscattering increases with 
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FIG. 5. (Color online) Simulated time traces of Eqs. Q5J) 
for the coupled setup shown in Fig. E] The upper (lower) 
trace shows the powers in SRL a (b). At t — 0.6 [is a trigger 
pulse is injected solely in SRL a [see Fig. |4j. Parameter 
values: 8 = 0.08, E t = 1.6 x 10" 4 , k c = 0.35fc, cj> c = 0.35tt. 
Conventions as in Fig. [T] 



the reflective coupling between the counterpropagating 
modes, which can be manipulated through the cavity de- 
sign. Due to the weak character of the coupling, excita- 
tory excursions still persist due to the similar phase space 
structure. We will, however, not project the trajectories 
on the respective (9, ip) planes of each SRL. The phase 
space structure of the solitary SRL in these planes does 
not provide a suitable reference since the projection does 
not capture the extra dimension of the coupled system. 

Coupling the two SRLs nevertheless has an effect on 
the dynamics of the individual lasers. Due to the sin- 
gle waveguide coupling SRL a and b are continuously in- 
jected in respectively the CCW and the CW mode, which 
stabilizes these respective modes in the SRLs [25]. On 
the other hand, the cavity asymmetry 8 favors the CCW 
mode in both SRLs. In the range of parameters we con- 
sider, the coupled system has a steady state s± that has 
both SRLs lasing in the CCW mode [see Fig. [I] , which is 
stable throughout the considered parameter range. We 
raised the value of 8 in this section to increase the sta- 
bility of this initial steady state. Furthermore, there is 
another stable steady state S2 in which SRL a and b are 
respectively lasing in the CCW and CW mode. This 
steady state does not exist for coupling phases 4> c close 
to (and 7r, due to the transformation in Eqs. ([6])). 

Fig. [5] shows that a pulse generated in SRL a — excited 
by an external optical trigger signal — can in its turn ex- 
cite a pulse in SRL b. Note that the pulse shape in SRL a 
is a little distorted. Small oscillations are superimposed 
on it and the damping to the resting state is increased. 
The superimposed oscillations are relaxation oscillations 
as explained in the previous section for high(er) ampli- 



tilde trigger pulses. The need for a slightly higher ampli- 
tude stems from the stabilizing effect of the continuous 
injection from SRL b, which also explains the increased 
damping rate toward the resting state. This change in 
damping rate is not very visible in SRL b [compare with 
Fig. [I] because the weakest mode of SRL a is injected 
in the CW (firing) mode of SRL b. The different power 
levels of the resting state in SRL a and b in Fig. [5] is 
an artifact of the single waveguide coupling, which pro- 
duces an additional asymmetry of the global system. The 
dominant mode of SRL b is coupled to SRL a, while 
the weakest mode of SRL a is coupled to SRL b. The 
dominant mode of SRL a and the weak mode of SRL b 
are not fed through. This gives rise to the asymmetric 
steady state power levels. If we would add an additional 
bus waveguide connecting the bottom of both rings, the 
dominant mode of SRL a and the weak mode of SRL b 
would also be fed through. No additional asymmetry is 
introduced in this case and the steady state power lev- 
els are symmetric. This double waveguide coupling also 
yields positive results with regard to pulse excitation in 
SRL b, but exhibits much more unstable behavior than 
the single waveguide coupling. For this reason, we will 
focus on the single waveguide coupling scheme of Fig. [4] 

The response of SRL a and b for various values of k c 
and c/> c are shown in Fig. [6] The color coding represents 
the number of times the line 9 — is crossed in the 
asymptotic (6, ip) phase plane, which corresponds to an 
equal power in both modes. Two subsequent crossings 
in SRL b can hence be regarded as an excited pulse. An 
odd number of crossings signifies that SRL b has switched 
from CCW to CW lasing, instead of returning to its rest- 
ing state and emitting a pulse. For the coupled system 
this indicates a switch to the steady state sa- This is 
the case for the pale blue (1 crossing), yellow (3 cross- 
ings) and part of the brown area (5 or more crossings) in 
Fig.|] 

Just as in the previous section, the optical phase plays 
an important role. Certain values of <p c will be favorable 
and others will be detrimental for transferring excited 
pulses. The role of <pc is however not the same as the 
role of the "phase difference" between the injected and 
the SRL fields in the previous section. When residing 
in their resting state SRL a and b will have a constant 
phase difference between their fields, which is a function 
of (j) c . Moreover, the phase of the "source" laser (SRL 
a) is no longer independent of the injected laser due to 
the bidirectional character of the coupling. Nevertheless, 
there is a region around 4> c = (and <j) c = 7r) for which no 
response can be elicited. Values of cj) c around ir/2 yield 
clear responses from SRL 6, i.e. the number of crossings 
is nonzero. The area with value 2, in which we excite 
a single pulse in SRL b, is relatively large. An example 
of a time trace in this area was shown in Fig. [5] and is 
indicated by a white circle in Fig. [6j The brown region 
at cj) c rs 0.6-7T indicates more complex behavior than sin- 
gle pulse excitations, an example of which is shown in 
Fig. [7] The corresponding location in Fig. [6] is indicated 
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FIG. 6. (Color online) The response of SRL a (top) and b 
(bottom) for different values of k c and cf> c . The colormap 
displays the number of times the line 8 = is crossed in the 
asymptotic (9, tp) phase plane. Number of crossings larger 
than 5 have been renormalized to 5 to improve the readability. 
The upper (lower) white line in the bottom figure bounds the 
region where the number of crossings in SRL a is > 5 (=4). 
The white circle (square) indicates the parameter values for 
the time trace in Fig. [5] (Fig. [?J|. Trigger pulse: E t = 1.6 x 
10 -4 . Parameter values: 8 = 0.08. 



by a white square. Note that both time traces have the 
same coupling amplitude but a different coupling phase. 
Instead of each firing a single pulse, the different value 
of 4> c allows SRL a and b to emit a series of spikes. Note 
that they are emitted in anti-phase by SRL a and b. Af- 
ter a while this oscillation can no longer be sustained and 
the SRLs return to their resting states. The range over 
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FIG. 7. (Color online) Simulated time traces of Eqs. (|5j 
for the coupled setup shown in Fig. B] The upper (lower) 
trace shows the powers in SRL a (b). At t = O.Qfis a trigger 
pulse is injected solely in SRL a [see Fig. [4]. Parameter 
values: E t = 1.6 x 1CT 4 , 8 = 0.08, k c = 0.35k, (j> c = 0.7n. 
Conventions as in Fig. [T] 



which both SRL a and b yield a complex response coin- 
cides with the existence of an unstable limit cycle. For 
increasing <p c it originates in a homoclinic bifurcation on 
the left side of the complex region (the side of lower (f> c ) , 
and is destroyed in a subcritical Hopf bifurcation on the 
right side of the complex region, resulting in the desta- 
bilization of S2. So different values of the coupling phase 
can have a profound influence on the dynamical behavior. 

Two coupled asymmetric SRLs can thus function as 
communicating neurons. We moreover observe that they 
do this in a single direction, i.e. the configuration of Fig. 
[4] enforces a preferential "propagation" direction for the 
excited pulses. We have shown that it is possible to excite 
a pulse in SRL b with a pulse from SRL a that was excited 
by an external trigger. However, when exciting a pulse 
in SRL b with an external trigger, we have not found 
any parameter values for which SRL b excites a pulse 
in SRL a. Note that in this case the "pulse" traveling 
from SRL b to SRL a is actually not a pulse, but rather 
a dip in the steady state power. The pulse can actually 
be interpreted as an increase in a destabilizing injection, 
while the dip is a decrease in a stabilizing injection [see 
Fig. [4] , explaining the difference in response sensitivity. 

V. CONCLUDING REMARKS 

In this work, we have theoretically investigated the 
possibility of triggering a pulse in an excitable (asymmet- 



ric) SRL by using an optically injected trigger pulse. We 
have used a standard rate equation model for the numeri- 
cal simulations and an asymptotic two-dimensional phase 
plane to interpret the results. This two-dimensional 
phase plane provides a well-defined threshold appearing 
as two branches of a stable saddle manifold that need 
to be crossed. Using this approach we have shown that 
the phase difference between the injected field and the 
SRL fields determines the direction of the perturbation in 
this phase plane. An unfavorable phase difference, which 
we were able to derive from the asymptotic phase plane 
structure, fails to excite a pulse in the SRL. Furthermore, 
a mechanism for exciting two or more consecutive pulses 
using a single trigger pulse was revealed. The appearance 
of such multiple response pulses is not related to period- 
doubling or homoclinic bifurcations as in [5J 110) , nor to 
noise-induced clustering as in [12] . They arise due to the 
multiple folding of the stable manifold, i.e. the threshold 
curve, around the resting state. 

We have also shown that two SRLs, coupled by a sin- 
gle bus waveguide, can excite each other and can thus 
function as communicating neurons. This type of neural 
network can be fully integrated on chip and does not suf- 
fer from the drawback of needing extra-cavity measures 
as other optical neurons do [2"HTT] . The typical pulse du- 
ration of 10 ns moreover yields a processing speed that is 
10 5 times larger than biological neurons. These advan- 
tages in size and speed, and the high degree of parallelism 
offer good perspectives as an artificial neural network. 

Although the asymptotic two-dimensional model used 
here is derived particularly for SRLs, the excitability 
mechanism presented in this work is general for a subclass 
of Z2-symmetric systems close to a Takcns-Bogdanov 
point providing the same sequence of bifurcations [13]. 
As an example of another optical system that will exhibit 
similar dynamics, we mention semiconductor disk lasers 
[2 6) as they essentially share the same circular symmetry 
as SRLs. 
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